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We calculate the cross section of a single inclusive gluon production in deep inelastic 
^ , scattering at very high energies in the saturation regime, where the parton densities inside 

hadrons and nuclei are large and the evolution of structure functions with energy is nonlinear. 
The expression we obtain for the inclusive gluon production cross section is generated by 
this nonlinear evolution. We analyze the rapidity distribution of the produced gluons as well 
as their transverse momentum spectrum given by the derived expression for the inclusive 
cross section. We propose an ansatz for the multiplicity distribution of gluons produced in 
nuclear collisions which includes the effects of nonlinear evolution in both colliding nuclei. 
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I. INTRODUCTION 



At very high energies corresponding to very small values of Bjorken x variable the density of partons in 
the hadronic and nuclear wave functions is believed to become very large reaching the saturation limit |^-^ . 
In the saturation regime the growth of partonic structure functions with energy slows down sufficiently to 
, unitarize the total hadronic cross sections. The gluonic fields in the saturated hadronic or nuclear wave 

' function are very strong [^^. A transition to the saturation region can be characterized by the saturation 

\ scale Ql{s), which is related to the typical two dimensional density of the partons' color charge in the 

infinite momentum frame of the hadronic or nuclear wave function [B-p|. The saturation scale Ql{s) is an 
increasing function of energy s and of the atomic number of the nucleus A ||7|-p^ . At high enough energies 
or for sufficiently large nuclei the saturation scale becomes much larger than ^^cd allowing for perturbative 
description of the scattering process at hand p.p^. The presence of intrinsic large momentum scale Qs 
justifies the use of perturbative QCD expansion even for such a traditionally non-perturbative observables 
as total hadronic cross sections. 

Recently there has been a lot of activity devoted to calculating hadronic and nuclear structure functions in 
the saturation regime. The original calculation of quark and gluon distribution functions including multiple 
rescatterings without QCD evolution in a large nucleus was performed in [^. The resulting Glauber-Mueller 
formula provided us with expressions for the partonic structure functions which reach saturation at small . 
McLerran and Venugopalan has argued in Q that the large density of gluons in the partonic wave functions 
at high energy allows one to approximate the gluon field of a large hadron or nucleus by a classical solution of 
the Yang-Mills equations. The resulting gluonic structure functions has been shown to be equivalent to the 
Glauber-Mueller approach |^,§|,|9| . An important problem which still remained was the inclusion of quantum 
QCD evolution in this quasi-classical expression for the structure functions. The problem was equivalent 
to resummation of the multiple BFKL pomeron [po[ exchanges. The evolution equation rcsumming leading 
logarithms of energy (a^ In s) and the multiple pomeron exchanges was written in using the dipole model 
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of g] and independently in |1C|] using the effective high energy lagrangian approach. The equation was 
written for the cross section of a quark-antiquark dipolc scattering on a target hadron or nucleus, which, in 
turn can yield us the F2 structure function of the target. The latter can be written as 
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(1) 



with ^''^(l, a) the wave function of a virtual photon in deep inelastic scattering (DIS) splitting in a 
qq with transverse separation r and the fraction of photon's longitudinal momentum carried by the quark 
a. ^''(zi, a) is a very well known function and could be found for instance in ||7|, p^ , p^ . The quantity 
N{r, 6, Y) has a meaning of a forward scattering amplitude of a dipole with transverse size r at impact 
parameter b with rapidity y on a target proton or nucleus normalized in such a way that the total cross 
section for the process is given by 



(fbN(r, b, Y). 
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The evolution equation for N closes only in the large- A^'c limit of QCD |jT^,0 and reads 



Nixo-i_,b,Y) = -7(xoi,6) exp 
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with the initial condition set by 7(3101 which is the propagator of a dipole of size Xj^^ the impact 
parameter b_ through the target nucleus or hadron. 7 was taken to be of Mueller-Glauber form in [Q: 
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where for a spherical nucleus ||3|,|19|,pi| 
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with p = A/[(4/3)7ri?^] the density of the atomic number A. The gluon distribution in Eq. (H) should be 
taken at the two-gluon order 
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with /i some infrared cutoff. The scale Qq^"^^'^^ has the meaning of the quasi-classical quark saturation scale 
generated by multiple rescatterings prior to the inclusion of evolution. Eqs. (|l|) and (|^) provide us with the 
F2 structure function and the total cross section of DIS on a nucleus including all multiple BFKL pomeron 
exchanges (fan diagrams). In spite of several attempts to solve Eq . (|3| ) analytically, which provided us with 
a well- understood high and low energy asymptotics for N fj Jl2| , |l l| , the exact analytical solution is still 
to be found. There exist several numerical solutions to Eq. ( |) demonstrating that at very high energies 
the amplitude goes to an independent of energy constant {N —> 1) thus unitarizing the total DIS cross 
sections [|l^,|l^,|l6|,0 . The numerical analyses also show that Eq. (||) does generate a momentum scale Qs 
which rapidly increases with energy. This scale justified the small coupling as expansion and helps avoid 
the problem of infrared instability of BFKL equation. An effort to calculate the NLO correction to Eq. (^) 
is currently under way p^ . 

Several other observables could be calculated in the framework of the saturation approach to hadronic and 
nuclear collisions. Of the exclusive observables diffractive (or, more precisely, elastic) cross section has been 
calculated in the quasi-classical limit in |23| . The evolution equation including multiple pomeron exchanges 
has been recently written for the cross section of single diffractive dissociation in pi] . 
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In this paper we are interested in inclusive particle production cross sections. In the saturation framework 
these cross sections have been extensively studied at the classical level in DIS, proton-nucleus (pA) and 
nucleus-nucleus (AA) collisions. The inclusive gluon production cross section for pA and AA in the weak 
classical field limit (lowest order in perturbation theory without multiple rescatterings) has been calculated 
in Ipsll , reproducing the result of Gunion and Bertsch [p6| . In the strong field limit including all multiple 
rescatterings but no QCD evolution the gluon production cross section was calculated for pA in [|l9| . This 
result has been recently reproduced in |27|-|29(|. The inclusive gluon production cross section for DIS was 
calculated in pH , with the resulting expression being slightly different from a straightforward generalization 
of the pA result of ||l9| . Finally an important problem for heavy ion physics is the calculation of the inclusive 
gluon production cross section in AA, which would provide us with the initial conditions for the possible 
formation of quark-gluon plasma. Numerical estimates of the related gluon multiplicities were performed in 
^0| , while an analytical ansatz has been proposed in |3l|| . 

The problem of inclusion of nonlinear evolution in the inclusive cross sections received much less attention 
in the literature. The case of heavy flavor production in DIS with nonlinear evolution has been solved in 
. The inclusive gluon production in DIS has been studied in |33|j3^ ] using the AGK cutting rules and 
fcy-factorization approach. 

In this paper we calculate the single inclusive gluon production cross section in DIS including the effects of 
multiple rescatterings and nonlinear evolution of Eq. (^. We begin in Sect. II by reviewing the derivation of 
the formula for the single inclusive gluon production cross section in the quasi-classical approximation given 
in ||2l[| . In the quasi-classical approximation the quantum evolution is not included since one is interested in 
resummation of multiple rescatterings . Each multiple rescattering of the produced gluon on a nucleon 
in the nucleus brings in a factor of a? A™ and the quasi-classical approximation can be defined as resumming 
powers of this parameter . The result for gluon production cross section is shown in Eq. ( p^ ) . 

We continue in Sect. Ill by including the effects of nonlinear dipole evolution from Eq. ^ in the expression 
for the cross section. The philosophy of our approach is similar to . We first construct a classical Glauber- 
Mueller-type expression for the inclusive cross section and then use it as our starting point for including 
dipole evolution. We are working in the rest frame of the target, which allows to consider the quantum 
evolution in energy as happening only in the wave function of the incoming qq pair [j7|pO|. In Sect. Ill A we 
analyze the evolution preceding the emission of the gluon that we measure in the final state. This evolution 
corresponds to emission of gluons with larger (harder) light cone component of momentum than the one 
carried by the gluon that we trigger. We show that this early evolution can only be linear (single pomeron 
exchange) leading to creation of the dipole in which our measured gluon is emitted. This conclusion is in 
agreement with the prediction of AGK cutting rules for inclusive cross section p3|-06|. In Sect. IIIB we 
proceed by analyzing the emissions of the gluons which are softer than the measured one. We demonstrate 
that the effect of these later emissions can be incorporated in the inclusive cross section by replacing the 
Glauber-Mueller exponents in Eq. (|l^) with the evolved forward amplitudes of a gluon (adjoint) dipole 
scattering on the target. The final expression for inclusive gluon production cross section in DIS is given by 
Eq. (1^) in Sect. IIIC which is the main result of this paper. 

We begin analyzing the cross section of Eq. ( pO[ ) in Sect. IV by observing that in the case of a very large 
nucleus corresponding to zero momentum transfer in each of the exchanged pomerons it can be rewritten 
in a factorized form as a convolution of two functions with Lipatov's effective vertex inserted in the middle 
(see Eq. (p9|)). In this form it almost agrees with the expression derived by Braun in |34] using the AGK 



3JI). To obtain the expression given in |Q 
amplitude of an adjoint dipole on the target 



rules and pomeron fan diagram approach (see Eq. (10) in 
from our Eq. (39) one has to replace in it the evolved forwarc 
nucleus Nq by the similar amplitude for the fundamental dipole N. While the difference seems minor in the 
weak field limit given by the linear evolution equations, it becomes much more profound in the transition to 
the saturation region where the quantities Nq and N obey different evolution equations while approaching 
the same high energy asymptotics. In Sect. V we study the transverse momentum spectrum and energy 
dependence of the obtained cross section. We observe that while at very large k± the cross section exhibits 
the usual behavior it softens to l/k'j_ in the small k± ~ Qs region similarly to the quasi-classical cross 
sections of |l9[ , ^ . The rapidity distribution of the produced gluons may have a maximum which position 
is determined by the values of the produced gluon's momentum k± and photon's virtuality Q (see Fig. |^). 

We conclude in Sect. VI by discussing the general principles of inclusion of nonlinear evolution in the cross 
sections for various inclusive processes. We present an ansatz for the multiplicity distribution of the produced 
gluons in nucleus-nucleus collisions (AA) including the effects of nonlinear evolution in both colliding nuclei. 
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II. INCLUSIVE CROSS SECTION IN THE QUASI-CLASSICAL APPROXIMATION 



In this Section we will review the derivation of the single inclusive gluon production cross section in DIS 
from pH ] in the quasi-classical approximation developed in [ p9| . The gluon production cross section in DIS 
can be rewritten as 



, I'A^qqGX 

dPk dy 
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(fxoi da $^ 



d^k dy 
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where da'^^^i / d'^ k dy{xQi) is the gluon production cross section for the scattering of a dipole of transverse size 
£01 = £0 ~ £1 '^^ target. We want to calculate this observable including all multiple rescatterings of the 
qq pair and the produced gluon on the nucleons in the target nucleus. In the quasi-classical approximation 
each rescattering happens via one or two gluon exchanges with the nucleon [ p^ . We assume that the initial 
quark-antiquark pair is moving in the "-f " light cone direction. The diagrams contributing to the gluon 
production cross section in the ~ light cone gauge are shown in Fig. ^ 




nucleus 



nucleons 





• • • 



FIG. 1. Gluon production in DIS in the quasi-classical approximation. The produced gluon may be emitted either 
off the quark or off the antiquark lines both in the amplitude and in the complex conjugate amplitude. Only one 
connection is shown. 
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Similar to [|T|,|ll|7|'|l,|l| the gluon emission can happen v a two possible scenarios: the incoming dipole 
may have the gluon fluctuation in its wave function by the time it hits the target or the gluon may be emitted 
after the dipole interacts with the target. Even with all the multiple rescatterings the interaction with the 
target is instantaneous compared to the typical emission time of the gluon . We thus may denote the 
interaction time r = x+ = 0. If the gluon emission time in the amplitude is ri and in the complex conjugate 
amplitude is T2 than the following classification of diagrams in Fig. |^ is possible. The graph in Fig. |i|A 
corresponds to the case when n < and T2 < 0, while the diagram in Fig. |i]b reflects the n < 0, T2 > 
case. A "mirror image" diagram should be added to Fig. |^ describing the ti > 0, r2 < case. The late 
emission ti > 0, T2 > scenario is represented in Fig. ^p. The produced gluon line can start off either quark 
or antiquark lines both in the amplitude and in the complex conjugate amplitude in Fig. 0. We are going to 
sum over all possible emissions while only one case is shown in Fig. |l|. 

To calculate the diagrams in Fig. |l| we will be working in transverse coordinate space with the intent 
to perform a Fourier transform into momentum space in the end. Thus the produced gluon has different 
transverse coordinates in the amplitude and in the complex conjugate amplitude {z^ and ^2)- The dipole 
consists of the quark at Xq and an antiquark at x^ . Due to real- virtual cancellations only the diagrams where 
the nucleons interact with the produced gluon survive in Fig. gA The square of the gluon's propagator 
could be easily calculated to give 

e-(£i-£2)'Qo./4 

with 



2^2 2 Sn'^asNc^/R'^ - h"^ ^, , ,2^ 
x^Ql, = ^ lifTx pxG{xM\x^). (9) 

The scale Q\g has the meaning of the gluon saturation scale in the quasi-classical (no evolution) approxima- 
tion and is different from the quark saturation scale of Eq. (^) by the Casimir operator. In Fig. |ip3 only the 
interactions with the gluon line and the antiquark are allowed. Performing the calculation along the lines of 
Appendix A in [p^ and adding the "mirror" contribution one obtains 

The minus sign in Eq. ( [lO| ) is due to the fact that the gluon is emitted after the interaction on one side of the 
cut in Fig. [^B. Finally the interactions with the nucleons shown in Fig. |^C would have canceled if we were 
considering proton-nucleus (or, more precisely, quark-nucleus) collisions as it was done in p9| . However, in 
the DIS case at hand the quark and antiquark originate from a virtual photon and thus have to be in the 
color singlet initial state on both sides of the cut. This condition was not imposed in pA where the 

color of the interacting quark was assumed to be randomized by the non-perturbative "intrinsic" quarks and 
gluons in the proton's wave function [ p7[ . Therefore, unlike the pA case moving an exchanged gluon line 
across the cut in DIS would modify the color factor of the diagram and thus real-virtual cancellation would 
not happen. A calculation of the interactions gives 

e-(£o-ii)^Qos/4 (XI) 

for the diagram in Fig. |l|C. 

To obtain the final answer we now have to combine the terms from Eqs. (||), (0) and (|ll|), multiply by 
the amplitude of gluon's emission in the original dipole while keeping in mind that the amplitudes in Fig. |l] 
are slightly different for different connections of the gluon to the qq pair. We then should Fourier transform 
the expression into momentum space. The final result for the single inclusive gluon production cross section 
of dipole-nucleus scattering in the quasi-classical approximation reads 



where h = (a;g -t-aij^)/2 is the dipolc's impact parameter. Together Eqs. (^ and ( [T^ ) give us the inclus 
gluon production cross section in DIS in the quasi-classical approximation as derived in [ pl| . 
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III. INCLUSION OF EVOLUTION EFFECTS 



We are now going to include quantum evolution into Eq. (|T^). We begin in Sect. IIIA by discussing 
the evolution preceding the emission of the measured gluon, and continue in Sect. IIIB by analyzing the 
subsequent evolution. The final expression is derived in Sect. IIIC. 

A. Emission of Harder Gluons 

Let us explore how the emission of gluons with the light cone "plus" component of the momentum much 
larger than that of the measured gluon modify the inclusive cross section. For simplicity we first consider 
an emission of a single extra gluon. The diagrams relevant for real emission of this extra gluon are shown in 
Fig. Ij. Since we are interested in the large- A'c dipole evolution all gluons are represented by double quark 
lines. Notations are explained in Fig. ^A. The incoming original dipole of transverse size Xjqi = Xq — SLi emits 
a (harder) gluon with transverse coordinate Then the measured gluon (#3) is emitted in one of the 
color dipoles formed by the emission of gluon ^2. Since we are interested in keeping this gluon's transverse 
momentum fixed in the final state its transverse coordinates are different on each side of the cut (xg and x^, ) . 
Emission in the lower dipole only is shown in Fig. |^. Emission in the upper dipole is completely analogous 
and could be obtained from Fig. || by switching gluons 2 and 3. 

ABC 








FIG. 2. Emission of a harder gluon in the dipole evolution. Gluons are denoted by double lines in the large Nc 
limit. The produced gluon is marked by a cross. 

The gluon lines in Fig. ^ can connect to either quark or antiquark lines in the dipoles off which they are 
emitted both in the amplitude and in the complex conjugate amplitude. This is shown in Fig. |2| by not 
connecting the gluon lines to any of the quark line specifically. For instance the gluon =f/=2 can be emitted off 
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the quark line and off the antiquark hne 1, which is demonstrated by drawing gluon jj=2 directly between 
those quark lines. This notation is the same as the one used in P,p^. The interactions with the target are 
not shown explicitly in Fig. |[ Instead we mark with a dashed vertical line in Fig. |^ the moment in light 
cone time r = when the multiple rescatterings of Fig. |l| occur (see Fig. |A). Of course the interactions 
may happen both in the amplitude and in the complex conjugate amplitude. The solid vertical line in Fig. ^ 
denotes the final state corresponding to r = cx). 

The measured gluon #3 is much softer than the gluon #2, fc3+ <g; fc2+, while their transverse momenta 
are not ordered. Thus the lifetime of the gluon #2 is much longer than that of gluon ^3. Therefore it seems 
natural that in all the diagrams considered in Fig. || the gluon #2 is emitted before gluon #3 on both sides 
of the cut. However the graphs in Fig. || do not cover all the diagrams which need to be considered. There 
is another set of diagrams, predominantly with the virtual emission of gluon #2, some of which are shown 
in Fig. ^ that could be important for gluon production at this order. We are going to first analyze the 
diagrams in Fig. ^. which would allow us to understand which ones of the diagrams omitted in Fig. |^ should 
be considered. 

The transverse momentum k and rapidity y of the measured gluon #3 are kept fixed, while the transverse 
momentum and rapidity of the other emitted gluon are integrated over. In the spirit of the leading logarithmic 
evolution integration over the rapidity of the gluon #2 is supposed to give us the factor of In 1 /x with x 
the Bjorken x variable, or, equivalently, a factor of total rapidity interval Y . This factor would make up for 
the suppression by the power of coupling constant Us that we have introduced by emitting that gluon. Now 
we are going to demonstrate that this enhancement does not happen in Figs. ^E-^. That implies that the 
diagrams in Figs. |E-|I do not give a leading logarithmic contribution and could be neglected as subleading. 




FIG. 3. One of tlie diagrams contributing to tlie class of graphs represented by Fig. 



Let us calculate the diagram in Fig. ^ which is one of the graphs represented by Fig. The diagrams 
in Fig. H are understood to be in the hght cone perturbation theory (LCPT) [Q, which has been used in 
the original construction of the QCD dipole model Therefore we will use the diagrammatic rules of 
LCPT in evaluating the graph Fig. ^ We are working in the frame where the quarks in the original dipole 
have large "+" components of their momenta while the nucleons in the nucleus have large "— " momentum 
components. Though only one rescattering is depicted in Fig. ^ our conclusions will be easy to generalize to 
include multiple Glauber rescatterings of Fig. |l]. The intermediate states are denoted by vertical dotted lines 
in Fig. 1^. In estimating the energy denominators of the intermediate states at t > one should remember 
that the "— " component of the nucleon's momentum also changes. According to the rules of LCPT the "— " 
component of the final state should be equal to the "— " component of the initial (incoming) state |^^. Thus 
one requires that for the diagram in Fig. |3| the following condition should be satisfied 



P- 



2k 



\2 



2+ 



2fc3 



(13) 



where we have used that for exchanged Coulomb gluons L = ||5|,|19|- Using Eq. (|13|) in evaluating, 
instance, the energy denominator of the rightmost intermediate state in Fig. ^ yields [ |38[ 



for 
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2fe2 + 



2fe3 + 



2fc3 + 



(14) 
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We have used the fact that fc3_|_ <C k2+ in the last approximation in Eq. (14). The rule for calculating the 
diagrams in Fig. |^ could be formulated as follows: the energy denominators for intermediate states with 
T > should consist of the sum of the "— " momenta components of all the intermediate gluons in the state 
minus the "— " momenta components of all the gluons in the final state. (Note that according to the rules 
of LCPT the "— " momenta components are not conserved at the vertices and all intermediate lines are on 
mass shell Q.) 

Using the outlined strategy one can show that the contribution of the diagram in Fig. ^ is proportional to 

dk2+ J] -'-^,2^7,2,. = / ^fc2+ -\.,-2.. ^ ~ 2^h2 , ~ (15) 

3+ 



3+ X2,\3 1^2 "-2+ •'«3+ \I^3J 112 "-2 



fe)'^.^2^2+ "''^3+ fe)^.^2^2+ fe)^Ai2^3+ k- 



with p+ a large momentum of one of the quarks in the original dipole (fc3+ ^ P+)- As can be seen from 
Eq. mm the dia gram in Fig. |^ does not contain any logarithms of energy. It should be compared to the 
contribution of, for instance, diagram in Fig. |2|A, which is proportional to 

rdk2+j-^ = T^ln^. (16) 



■3 + 



k2+ kz+ k3+ k3+ 



This diagram is enhanced by extra logarithm of energy (or, equivalently, an extra factor of rapidity) and 
should be included in our leading logarithmic analysis. In the same approximation the diagram in Fig. ^ is 
subleading since it does not have a logarithm of energy in it. Similar calculations can be carried out for all the 
diagrams in Fig. |^ showing that they are subleading and should be neglected. Finally the same conclusion 
will be reached if one analyzes all diagrams in Figs. |1E-||i : they do not bring in an extra logarithm of 
energy and therefore must be neglected. 

Here we would like to remind the reader about the approximation we are using. It is the same approxi- 
mation as outlined in for the calculation of total DIS cross section. The (nonlinear) quantum evolution 
resums the powers of agNcY and is a function of this parameter (1^ is the rapidity variable.) At the same 
time we are interested in resumming multiple rescatterings, which also bring in extra powers of ag P, p^|j3l| ]. 
The multiple rescatterings are limited to two gluon exchanges with each nucleon. Thus the actual parameter 
which is being resummed by multiple rescatterings is a^A^/"^ [p yis|j3l[ |. Therefore the combination of leading 
logarithmic large- iVc evolution and Glauber-Mueller multiple rescatterings resums all powers of both agNcY 
and c?gA^I'^ in the cross sections. From this point of view the diagrams in Figs. ^E-^ are suppressed by an 
extra power of a^iVc not enhanced by an extrapower of Y . 

We have proven that of all graphs in Fig. || only diagrams in Figs. HA-^D contribute. Based on this 
information one may conjecture that the evolution preceding the emission of the measured gluon consists 
of gluons emitted before interaction both in the amplitude and in the complex conjugate amplitude. To 
complete the proof of this statement we have to consider a whole class of diagrams where the line of gluon 
#2 is shorter or of the same length as the line of gluon #3. These diagrams include virtual diagrams where 
gluon #2 is not present in the final state. Many of these diagrams are also subleading, similar to Figs. ||E-||[. 
This can be seen by performing the analysis outlined above and using the rule presented in the Appendix of 
for calculation of virtual contributions. At the end we are left with the diagrams shown in Fig. ^, each of 
which gives a leading logarithmic contribution. Not all of the diagrams in Fig. ^ are symmetric with respect 
to horizontal (left-right) mirror reflections. Therefore one should add mirror images to diagrams A-D, F-P 
which will be denoted by primes (e.g. A', B', etc.). We note that due to momentum conservation the gluon 
^2 can be emitted only off the quark lines and 1 in the original dipole. Therefore moving part or all of 
gluon #2 across the cut does not change the transverse coordinate (and, therefore, transverse momentum) 
structure of the diagrams. Using the cancellation of final state interactions demonstrated in ||] we observe 
that in Fig. | 

B + C = B' + C' = (17a) 
D + E + D' = (17b) 
H + I ^ H' + r = (17c) 
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J + K = J' + K' = 
L + M + N = L' + M' + N' = 



P + Q + P' = 0. 

We are left only with the diagrams A, F, G, O and their mirror images A', F', G', O' in Fig. ^. 
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(17d) 
(17e) 
(17f) 



FIG. 4. Another set of graphs contributing to emission of a harder giuon in the dipole evolution. 
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Combining the results of the analyses of the diagrams in Figs. ^ and ^ we conclude that only the diagrams 
where gluon #2 is either emitted or both emitted and absorbed before the interaction at r < survive. The 
non-vanishing diagrams are graphs A-D in Fig. || and A, F, G, O, A', F', G', O' in Fig. ^. One can easily 
check that these diagrams add up to give us one rung of the linear dipole (or, equivalently, BFKL) evolution. 
The diagrams A-D in Fig. |^ provide us with the real part of the dipole kernel, splitting the original dipole 
01 in two, in one of which we continue the evolution by emitting the measured gluon ^3. Diagrams A, F, 
G, O, A', F', G', O' of Fig. ^ yield us with virtual corrections to the same process. 

Now we are in a position to generalize our conclusions to the case of many harder gluons in the virtual 
photon's wave function. The evolution preceding the emission of the measured gluon is the usual early-time 
dipole evolution leading to creation at early times (before the interaction, r < 0) of the dipole in which the 
measured gluon is emitted. The effect of this evolution is to modify the probability of finding this dipole in 
the virtual photon's wave function. 

This preceding evolution can only be linear. This can be understood by analyzing diagrams A-D in Fig. ^. 
Real emissions in the dipole model may lead to branching of one dipole evolution into two simultaneous 
evolutions, being thus equivalent to triple pomeron vertices in the traditional language |f7|,p|, p^|j39| ] . For 
instance, in Figs. ||A-pI) there could in principle be some subsequent evolution in the dipole 02 leading to 
creation of many dipoles which would interact with the target. The evolution may also lead to creation of 
dipoles at late times (r > 0). However all these emissions and interactions would not modify the momentum 
of the gluon that we measure, since they happen in a different dipole isolated from ours in the large-A^c 
limit. Therefore one can show that all these interactions would cancel via real-virtual cancellations: first 
one can show the cancellation of Coulomb gluon exchanges similarly to some of the cancellations in Sect. II 
and then the cancellation of evolution in dipole 02 follows due to probability conservation (see ^,^). This 
can be done for any graph where the evolution branches into two with the emitted gluon being produced by 
one of the subsequent evolutions. We therefore conclude that the early evolution is linear. This conclusion 
is similar to what one would obtain applying AGK cutting rules to the process ||3^-|3^. We will return to 
this similarity in the next Section. 

To include the preceding evolution into the gluon production cross section of Eq. (0) one should thus 
substitute 

/ xoi da ^''"^'^^Xf^^, a) ^ [ d^rda^''''^'i^r,a) d'^B n{r,Xf^^,B- b,Y - y) (18) 
J J 2nxQ-^ 

where n{r,XQi,b,Y) is the number density of dipoles of size xqi at an impact parameter b where the mo- 
mentum fraction of the softest of the two gluons involved in making up the dipole is greater than . The 
quantity n(r,XQi,b,y) was defined in 1^ and was shown to satisfy the linear evolution equation equivalent 
to BFKL. The solution of that equation gives Q 

J d'Bn{r,x,,,B,Y) = J ^ e-^^^ 1^^)'', (19) 

where xW is the eigenvalue of the BFKL kernel ||2^ defined as 

X(A) = 2V(1)-^(1-A)-7^(A) (20) 

and 



(21) 

Integration in Eq. (^9|) runs along a straight line parallel to imaginary axis to the right of all the singularities 
of the integrand. In Eq. (|l^) Y is the total rapidity interval of the DIS process while the measured emitted 
gluon is located at rapidity y. 



B. Emission of Softer Gluons 



In the previous Subsection we have understood how to include the evolution of the harder gluons in the 
virtual photon's wave function. Here we will address the question of how to include the evolution of gluons 
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with the hght cone component of their momenta being softer than the measured one. For instance, in Fig. g 
after being emitted the measured gluon spUts the dipole 02 into two off-forward dipoles. The dipoles are 
off-forward because the transverse coordinate of the measured gluon is not the same in the amphtude and in 
the complex conjugate amplitude. Nevertheless there could still be dipole-like evolution in each of these two 
dipoles. More gluons with light cone components of their momenta much smaller than k^^ can be produced 
and these gluons can interact with the target as well. 

We claim that in order to include the effects of this quantum evolution (in the leading logarithmic approx- 
imation) into Eq. (^2|) one has to substitute 

l_e-i'Qo.®/4 NG{x,Ly) (22) 

for all the Glauber exponents in it. The quantity Ncix^b^y) is the forward amplitude of a gluon dipole 
scattering on the target. This quantity is similar to N{x,b,y) except the initial state for Nq consists of 
a pair of gluons in the color singlet state instead of a quark-antiquark pair. The normalization of Nq is 
analogous to Eq. (^. In the large- iVc limit an adjoint (gluon) dipole can be decomposed into two fundamental 
(quark) dipoles. Therefore the scattering amplitude of a single adjoint dipole on a target nucleus in the large- 
Nc approximation is equivalent to the scattering of two fundamental dipoles on the same target. One can 
thus conclude that 

NGix,b,y) = 2N{x,b,y)-N^{x,b,y). (23) 
The evolution equation for NG{x,b, y) can be obtained by inverting Eq. ( p3|) to express N in terms of Nq 

Nix, 6, y) - 1 - ^Jl-N^{x,b,y) (24) 

and by substituting Eq. (|^ into Eq. (^). 

One can see right away that without evolution {y ~ 0) Eq. ( p2| ) turns into equality (see, for instance, Q) 

^^G(ai,6,0) = 1-6-2'^'=®/^ (25) 

Our goal now is to prove that the substitution of Eq. (^) does correctly incorporate the effects of subsequent 
evolution in the gluon production cross section. As we now know the preceding quantum evolution happens 
only at r < and does not interfere with the subsequent one. Similarly to Sect. II one should consider four 
different cases, corresponding to different emission times of the measured gluon in the amplitude (ti) and 
in the complex conjugate amplitude (T2). The four cases are: (1) ti < 0, r2 < 0; (2) ri < 0, r2 > 0; (3) ri > 
0,T2 < 0; (4)ri > 0,T2 > 0. 

Let us first consider case (1). In the quasi-classical case of Eq. (|^) it gave us the exponent in 
Eq. (||). According to Eq. ( p2| ) we now have to prove that evolution replaces that exponent with 
1 — Ng (zj^ — + z_2),y) ■ That means that the gluons in the subsequent evolution connect only to the 

emitted gluon line mimicking the scattering amplitude of a gluon dipole of size z-^ — Z2 on a nucleus. 




FIG. 5. Diagrams including one softer gluon produced after the emission of the measured gluon. 
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To prove that the gluons in the subsequent evolution connect only to the gluon line of the measured gluon 
let us again consider a simple case when we add one more gluon emission to the quasi-classical diagrams of 
Fig. |l|. Different from the previous Subsection the extra gluon's light cone momentum should be much softer 
that the measured gluon's one. To prove the cancellation of all the diagrams where the extra soft gluon 
interacts with the quark lines it is sufficient to consider graphs presented in Fig. |^. There after emitting the 
measured gluon #3 we include and extra emission of a softer gluon #2 (^2+ <SC A;3+). All the other relevant 
diagrams could be obtained by vertical and horizontal reflections of the graphs in Fig. |^. In Fig. ^ we kept 
only the diagrams where the emission of gluon #2 is enhanced by a logarithm of energy leaving out all the 
subleading ones similar to how it was done in obtaining Fig. ^. The diagrams where the soft gluon (#2) 
is both emitted and absorbed by quark lines are easily canceled by real-virtual cancellations and 
need not be considered in much detail. In the diagrams A and F in Fig. |^ we imply summation over both 
possible time orderings of gluon ^2. As in Fig. ||the interaction with the target is not shown explicitly and 
is denoted by dashed vertical lines. 

Let us consider diagrams A and D in Fig. 0. Emissions of gluons #2 and #3 bring in the same transverse 
coordinate dependence in all graphs in Fig.^. The only possible difference between A and D could be in 
the interactions with the target. In diagram A the Coulomb gluons could be exchanged only between gluon 
#3 and the target, while it appears that in diagram D the target could interact with gluon #2 as well. 
(Interactions of the target with the quark and antiquark lines cancel via real-virtual cancellations similar 
to Fig. I^A.) However, since the transverse momentum of gluon #2 is being integrated over the transverse 
coordinates of that gluon are equal on both sides of the cut. That way the Coulomb exchanges between the 
target and this gluon also cancel through real-virtual cancellations. Therefore both graphs A and D include 
the same interactions with the target giving identical absolute contributions to the cross section. The only 
difference is that in graph D the gluon #2 is a real gluon while in graph A the gluon ^2 is purely virtual 
and is emitted and absorbed at t < 0. Thus the contribution of the diagram in Fig. ppV comes with the 
negative sign with respect to the diagram in Fig. ||d. Therefore they cancel each other: 

A + D ^ 0. (26) 

In considering diagrams B and C in Fig. | we note that the interactions with the target are manifestly the 
same in both of them. The only difference between graphs B and C is that the gluon has been moved across 
the cut in C. Using the cancellation of final state interaction rules outlined in ||^ we argue that these two 
diagrams cancel each other 

B + C = 0. (27) 

Similarly one can show that 

F + G = 0. (28) 

Finally a simple analysis can show that in graphs E and H in Fig. || the interactions with the target involve 
only gluons and are identical. In both cases the interactions are equivalent to scattering of a dipole 32 and a 
dipole 3'2 on the target. Thus the diagrams E and H give the same absolute contributions. However, their 
signs are different since E is virtual while H is real. Therefore they also cancel 

E + H ^ 0. (29) 

We have proved that only interactions with the gluon #3 survive in the symmetric case (1). The result could 
be easily generalized to any number of softer gluons. 

Lastly we have to show that the softer gluons that connect to 3 and 3' add up to give a scattering 
amplitude of a gluon dipole 33' on the target nucleus Nci^ ~ 3'). To prove this one has to carefully analyze 
the diagrams with one extra gluon, similar to Fig. || only with the gluon #2 connecting only to gluon #3 
(#3') and then generalize the result to any number of soft gluons. This is what has been checked by the 
authors. Alternatively we can use the duality property of the amplitude which was used in |Q. We may just 
argue that the amplitude would remain the same if we mirror-reflect the gluon #3' in the complex conjugate 
amplitude into the amplitude. Then the gluon production graph would become manifestly equivalent to the 
scattering of a 33' dipole on the nucleus justifying the substitution of Eq. (p2[) in case (1) as desired. 
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The proofs of substitution of Eq. (|2^ ) for cases (2), (3) and (4) is analogous to the one outUned above. In 
cases (2) and (3) one first has to show that emission of a softer gluon is only possible off the gluon and the 
quark involved in the appropriate exponents in Eq. (|lo|). That means all interactions with one of the quarks 
should cancel. Which one of the quarks becomes the "spectator" depends on the way the gluon #3 was 
emitted in the amplitude and in the complex conjugate amplitude, similar to Sect. II. After that by moving 
(reflecting) the remaining (interacting) quark across the cut one can show that, since in the large- A^c limit 
two quarks with the same coordinates are identical to a gluon, the interaction is identical to a scattering of 
an adjoint dipole on the target. The adjoint dipole would be composed of the gluon #3 and the interacting 
quark, thus justifying the substitution of Eq. (p2[). In case (4) we explicitly have two fundamental dipoles 
01 on both sides of the cut developing the evolution and interacting with the target. Using Eq. ( |23| ) we can 
once again prove the substitution of Eq. ( |2^ ) . 



C. Expression for the Inclusive Cross Section 



Now that we have understood above the effects of evolution on the quasi-classical expression for the 
inclusive gluon production cross section we can combine the results to write down the answer including all 
evolution effects. The evolution preceding the emission of measured gluon can be included by the substitution 
shown in Eq. (|l8|). The evolution is linear, similar to what one would obtain from AGK cutting rules ||3^-|3^. 
This is not the first time AGK cutting rules are shown to work for an observable in dipole evolution. They 
were also satisfied by the equation for the diffractive structure function found in |2j] . 

The evolution following the emission of the measured gluon is nonlinear and can be included by making the 
substitution of Eq. ( ^ ) . It is quite surprising that such a complicated evolution effect can be incorporated 
by such a rather compact rule. 

Combining Eqs. (^ and ( |l^ ) with the prescriptions of Eqs. (|lj) and (^2|) we obtain 



(Pk dy 



1 

2^ 



d^rda^^'-^''''{r,a) d^B n{r,Xj^^,B-b,Y^y) 



d^xoi agCp 1 



■d^b 



X d zid Z2& 



+ ^^G ( ^2 - o U2 + 2/ ) - ^G ( a - o Ui + I2). y ]-'^G\x^- Xj, -{Xi + Xj), y 



where the produced measured gluon has rapidity y, the total rapidity interval in Y and 



(30) 



(31) 



Eq. ( pO| ) gives us the single inclusive gluon production cross section in DIS on a hadron or nucleus including 
the effects of nonlinear leading logarithmic evolution and multiple rescatterings. This is the main result of 
the paper. 



IV. FACTORIZED FORM OF THE INCLUSIVE CROSS SECTION 

We now continue by analyzing Eq. (^0|) in the limit of a very large target nucleus. In that case the 
momentum transfer to the nucleus is cut off by inverse nuclear radius and is very small. We can therefore 
take the scattering amplitudes Nq in Eq. ( ^0| ) at t = 0, which in coordinate space is equivalent to neglecting 
the small compared to the nuclear radius shifts in impact parameter dependence. Similar to how it was done 
in 0,|l^ we assume that 

NGix,b±y/2,Y) ^ NGix,b,Y), (32) 
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where y represents any of the impact parameter shifts in Eq. (|3C|). Integrating over or z_2 depending on 
the argument of Nq in the term involved we reduce Eq. (^0|) to 



(Pk dy 



27r2 



fd^rda $^*^«« (r , a) d^B n(r , Xo^,B-b,Y-y) — 



■ik-z 



In- 



1 



\z-XiA^ \z-x,AA 



271 ^ 5 (z — X,. 



dHNG{z,b,y) (33) 



with x^j — x^ — Xj and A some infrared cutoff put in to regulate the integrals |^ . Performing summation 
over i and j, employing Eq. ( |l9|) and integrating over xoi we obtain 



^^j'A-,qqGX I 



27r2 



d^k dy 



/ TTl \zJ TT 



'r\2A a,Cp 1 



27r 



2i k 



1 

2A2 



1 2 



(34) 



Eq. ( p4| ) can be rewritten as 
d^fc 0??/ 



^ |d2.da$-*--(r,a)^^i-^ ld'zNa{z,y) 



X v: 



7ri Vz/ A2 



(35) 



with 



NgU.v) = / d'foA^G(^,fo,y)- 



(36) 



Integrating by parts yields 

^^-r'A^qqOX I 



d^k dy 



— / d^da^'' """'(i:,a 



Cf 1 2 



7r2 27r /c 



/ d?ze-^^-^ \^lNG{z,y) 



X (^d'Bnir.z^B^Y -y) 



Eq. (^) can be rewritten with the help of Lipatov's effective vertex 



AT 4_2 



(37) 



(38) 



as 



^^-r'A^qqCX I 



d^k dy 7r(27r)i 



d^rda $'^'^««(r, a) 



X / d2zArg(z,y)4(z)ij (^^d2Bn(r,z,5,y-y) 



(39) 
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Eq. ( p9| ) presents our cross section in some sort of fcr-factorized form, similar to the inclusive cross sections 
derived in |^^. It is almost identical with Eq. (14) in jsj], which was obtained by applying AGK rules to 
BFKL pomeron fan diagrams in DIS case. The only difference between our Eq. ( ^9|) and Eq. (14) in is 
that the latter uses N{z,b, y) (the function $(y, r, b) in the notation of Q) instead of our Nc{z,h, y). That 
is in order to obtain Eq. (14) of from our Eq. ( |39| ) one has to substitute Nq with TV in it. In the weak 
field limit when the evolution is linear both N and Nq obey the same BFKL evolution equation, though 
with different initial conditions given by Eqs. (^ and ( psf ) correspondingly. In the saturation region the 
evolutions equations for the two quantities in question become different. N{z,b,y) obeys Eq. (^, while the 
equation for Nq (z_, b, y) can be obtained by substituting Eq. (|24 ) into Eq. (||) . The resulting equation would 
involve square roots and would be quite different from Eq. ^). However the high energy asymptotics of 
the two objects is still the same: both N{^,b, y) and NG{z,b, y) asymptotically approach 1 at high energies, 
when y oo. 

The form of the cross section in Eq. (^9|) suggests that in a certain gauge or in some gauge invariant way 
it could be written in a factorized form involving two unintegrated gluon distributions merged by effective 
Lipatov vertex. The usual form of the factorized inclusive cross section is 

da 2a, f. /i(a,g^)/2(6Jfc-gn 

d q =^ ^— (40) 



d^kdy Cpk^ J f(fc-# 

with ^1 and ^2 the values of Bjorken x variable for the gluons in each of the colliding particles. /i('^i,(?^) 
and f2{£,2, \k — gp) are the unintegrated gluon distributions 



d xGi (£i , , 



Eq. ( pO| ) can be written in the coordinate space as 

da 27r2 



d^k dy NcCf 

with 



d'zMii,z)Lkiz)Mi2,z) (42) 



(27r)2 ^ qi 



Mi^,^) = / 7wJ:^e-'-^^^^- (43) 



Comparing Eq. ( |4^ ) with the gluon production cross section in dipole-nucleus scattering which follows from 
Eq. (§) 

^^qqA^qqGX ^ /" ,o , , . . - _ 1 / 1 



d^kdy ^,,)Jd^^^oi^^y)M,)^^{^d^Bnir_,^,B,Y^y)) (44) 

and assuming that the latter was generated through the same factorization mechanism we identify 

NG{z,y ^ ^ ^^^^cl^itz). (45) 

Employing Eqs. ( ^l| ) and (^) we may rewrite Eq. (^5|) as 

which corresponds to the definition of gluon distribution used in |l^,^^. However Eq. ( ^6|) is not satisfied 
in the quasi-classical limit of no evolution. There A^g(^, 0) is given by Eq. ( |25| ) integrated over impact 
parameter 

NgU, y^O) ^ S^{i- e-^- ^-Z^) , (47) 
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while the gluon distribution including all multiple rescatterings have been calculated in |]l9| to give 

r^,^^,,.,d^G^|^ = 2 /^2^TrM^^(0).A^^(z)) = (i-e-^^Q-/4), (48) 

where A^^ (z) is the non-Abelian Weiszacker- Williams field of a nucleus |^ and we assumed that the nucleus 
is cylindrical with the cross sectional area S± = ttR^. Expanding both Eq. (^) and Eq. (^) to the lowest 
order in Qn^ corresponding to two gluon exchange we can see that Eq. ( p6| ) can be easily satisfied. However 
full Eqs. (|7|) and ( psf ) when inserted into Eq. ( ^q ) do not satisfy it. Therefore Eq. (|6|) seems to work for the 
leading twist two gluon exchange approximation but appears to fail once we include multiple rescatterings 
in it. 

The failure of Eq. (^6|) in the quasi-classical limit makes the physical meaning of factorization of Eq. (^) 
quite obscure. The factorized form of Eq. implies convolution of two unintegrated gluon distributions 
with Lipatov vertex as shown in Eq. (^0|). It appears, however, that we can not identify one of the convoluted 
functions {Ng) with the unintegrated gluon distribution in disagreement with the factorization hypothesis 
of Eq. (^) . This observation by itself would be quite natural indicating that higher twists in the form of 
multiple rescatterings modify the relationship between Nq and the unintegrated gluon distribution. On the 
other hand the fact that fcy-factorization is still preserved seems somewhat unexpected. The reason why we 
are able to write down the inclusive cross section of Eq. (|3^) in the factorized form of Eq. needs to be 
better clarified which will be done elsewhere. 



V. PROPERTIES OF THE INCLUSIVE CROSS SECTION 

It was shown in that asymptotic solutions to the nonlinear evolution equation bear most of quantitative 
features of the general one and are very convenient for discussion of various high energy processes. To study 
asymptotic behavior of the inclusive gluon production it is convenient to integrate over directions of vector 
z in Eq. (^) and use Fourier transform of the gluon scattering amplitude Nail, v) defined as 

/•OO 

Naiz.y) ^ z^ dl I Jo{lz)NG{l,y) (49) 
Jo 

with / = |/|, obtaining 




W Q % 

FIG. 6. Various asymptotic kinematical regions discussed in the text. Solid line is the saturation scale Qs{y)- Bold 
lines are the boundaries of different regions. 
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There are four interesting asymptotic regions of Eq. (^o|). They are shown in Fig. ^. In the first and 
second regions the produced gluon momentum k is the largest external momentum scale of the process. The 
first region in Fig. || corresponds to the cases when k ^ Qs{y) ^ Q and k » Q :s> Qs{y), while in the 
second region Qs[y)'^ k^ Q and Qs{y) ^ Q ^ k. We can find asymptotic behavior of the inclusive cross 
section by expanding function xW near its simple pole at A = and evaluating the Mellin transform in A in 
the saddle point approximation. This corresponds to the double logarithmic approximation to the evolution 
equation which means summation of a" In" k {Y — y)" ~ 1 terms. The third and fourth regions in Fig. ^ 
depict the cases when Q S> fc 3> and Q'^ Qs{y)^ k correspondingly. In the third and fourth regions 

asymptotic behavior of the inclusive cross section can be found by repeating this procedure near the simple 
pole of x(-^) at A = 1 which again corresponds to summation of a" In" k {Y — y)" ~ 1 terms. These are two 
asymptotic double logarithmic regions relevant to the linear evolution before the gluon has been emitted. 
Generally, after its emission the evolution is nonlinear. So, there is a kinematical domain where the parton 
density is large (regions 2 and 4 in Fig. ^). In this domain a color dipole evolves into two dipoles, one of 
which is of the size 1/Qs{y)- And there is a kinematical domain where the parton density is small (regions 
1 and 3) and the evolution is still linear. 

Integration over z in Eq. ( po| ) yields a combination of generalized hypergeometric functions which can be 
expanded near A = and A = 1 to give 



dz /r\2>' 



(-) z^Mlz) 

Q Z \Z) 



2kz k'^z'^ 

— Mkz) - -^Mkz) + 2Jo(fcz) 



4(fcr)^- 



A ^ 0, A: > 



-, A 0, /c < 



(51) 



/c^(l-A)' A ^ i, ^ i, 



A 1, A; < Z. 



Each of the limits in Eq. (^ij) corresponds to the appropriate region in Fig. ^. It is now straightforward to 
evaluate the Mellin transform in the saddle point approximation, which corresponds to the double logarithmic 
approximation of the scattering amplitude. Using Eq. ( pl| ) in Eq. ( |50| ) we obtain in the region with k ^ I 
and A = 



cPkdy 



(2^)3 TT Ai^ 



g2^a,(y-a)ln(fcV2) 

{as\n{k^r^){Y - y))^^'^ Jo 



dllNaily), (52) 



We used the well-known expansion x(-^) = +0{\). Scattering amplitude Nq{1, y) must be normalized to 
give the correct Fourier transformed initial condition Eq. (^7|). We have (neglecting logarithmic dependence 
of Qs on I) 



91 



dX / I 
27ri \Qo 



2A 



,-2A-l 



goy^-^ r(i-A) 1 / £_ 

Q^s) l-A 2 \ 'Ql 



(53) 



where r(0, z) is the incomplete gamma function of zcroth order. In the region 1, where I ^ Qg, evolution is 
linear and the scattering amplitude is small. In the linear region 



Ncily) 



QlgirR^ f dv 



2P 



I 



2-Ki \Qo 



'\^y ^ Ql^f V^^l^'y'[\ 2^..yMiyQl)^ ln(//go) » a.y, (54) 



2nP 



where we expanded xi^) near v = and assumed the nucleus to be a cylinder with the cross sectional area 
S± = ttR^. Qq is an initial scale for linear evolution. Eq. ( |5^ gives gluon scattering amplitude in double 
logarithmic approximation with the expansion parameter as y\n{l/Qo) ~ 1. Upon substitution of Eq. ( p^ ) 
into Eq. (^2|) and integration over I one arrives at the asymptotic expression for the inclusive cross section 
in the region 1 
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(Pk dy 



region i 



(27r)4 2fc^ J ay"[y(y-y)]i/4(lnfc2/Q2ji/4ini/4(p^2) 

(55) 



In the region 2 with / <C saturation sets in and the total scattering amphtude N{x, y) is close to unity. 
Thus, using Eq. ( p3| ) we conclude that NQ{x,y) is close to unity as well. In the momentum representation 
this is equivalent to the following asymptotic behavior in the saturation regime |l2| 

NG{l,y) = HQs{y)/l), l<.Qs{v). (56) 

Using Eq. ( |56| ) together with the second line of Eq. ( ^ in Eq. (|5|) yields 



d^kdy 



region 2 



(2^)4^ fc2 y V-' ^(c,4r_y)ln(fc2^2))l/4- 



Consider now A = 1 pole of function x(A) in Eq. (^o|). In the vicinity of this point x(A) = (1 — A)~^ + 
0(1 — A). Substituting the third line of Eq. (^ij) into Eq. ( |5l| ) and repeating the above procedure we end up 
with following asymptotics in the third kinematical region. In the region 3 (Q S> A; S> Qs) 



7r/?2 n r , 2^a,(l'-y)ln(l/fc2r2)+2Va.yln(fcVQL) 



regions (2^)^ F J ' aJ^[yiY - y)Y/^\n{kyQlJ\n{l/k^r^))^^^' 

(58) 



(i2fc fijy 



Inserting the last line of Eq. (|5l|) into Eq. (^0|) and using Eq. (|5^) together with the expression for Nq 
generated by linear evolution gives for the region 4 (fc <C Qs <C Q) 



d^^'A^qqGX 



d^kdy 



region 4 



47ri?2a,02() , ^2Va.(y^,)i„(i/Q?(.).^) 

-—7 , ' Vtt / d'^rda^^ '^'^(r,a)r-' — , (59) 



where the integration over I has been done with logarithmic accuracy. 

Note that the spectrum of produced gluons falls off as ^ k~'^ in region 1, which is a well-known result of 
the perturbation theory. One factor of fc~2 arises from the Lipatov's effective vertex, while another one is 
given by the perturbative behavior of the scattering amplitude V'^Ng which scales as ^ fc~2 at large fc_L. 
When the typical momentum / of the late evolution is less than the saturation scale Qs{y), the anomalous 
dimension of the gluon structure function is close to unity. Thus the scattering amplitude depends only 
logarithmically on momentum. As a result, in the region 2 the gluon spectrum softens to fc~2 behavior. 
In the region 3 the late evolution is linear like in region 1. However an additional factor of fc2 stems from 
different double logarithmic asymptotics at A = 1. It is remarkable that unlike the regions 1 and 2, the 
spectrum in region 4 drops in the same manner as in the region 3 despite the fact that the late evolution is 
nonlinear. This happens because the typical value of I in the rightmost integral in Eq. ( |5^ is / ~ Qs- Due 
to this fact Eq. (^) is enhanced by a factor of rather than fc2 Qf course, one could not expect that 
the spectrum will depend on produced gluon's momentum only logarithmically in the presence of a hard 
projectile (virtual photon) which wave function has not reached saturation yet. 

Derived asymptotic formulae Eq. (|55|), Eq. (p7|), Eq. (^8|) and Eq. ( |59| ) make it possible to understand 
some important features of the spectrum energy dependence. Equating the y derivative of the expression 
in the exponent of Eq. ( |55| ) (or Eq. (|5^)) to zero one finds the position of the inclusive spectrum maximum 

at some fixed k in the double logarithmic approximation in which the pre-exponcntial factor is a slowly 
varying function. It reads 

Vo = y , , — — > -Y, (60) 

ln^(fe/Qo,) 

where -I- and — correspond to the regions 1 and 3. By analogy we find that in the region 2 spectrum is 
a monotonically decreasing function of y and exhibits no maximum. In the region 4 the spectrum has no 
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maximum in y as well. The following equation defines a surface in the space of parameters where spectrum 
is constant: 



AoisY + In 



;(0)r 



= 16a,yin 



Q2(0)r 



(61) 



Although Eq. ( |6l|) is beyond the validity of the double logarithmic approximation employed in this section, 
it shows the general feature of the spectrum y-dependence in region 4. Namely, when Q ~ Qs the spectrum 
is a step-like function. In Fig. ^ we show a qualitative picture of the energy behavior of the inclusive cross 
section given by Eq. ( |30| ) in different kinematical regions. 
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FIG. 7. Qualitative energy dependence of the inclusive gluon production cross section. All lines are scaled by 
different numerical factors to underline important features of the spectrum. In all plots we show the curves for two 
different values of k and/or Q (fci < k2, Qi < Q2)- 

Note, that all four plots in Fig. ^ show different qualitative behavior of the inclusive cross section as 
the rapidity varies. By plotting the experimental data on the energy dependence of the spectrum one can 
distinguish between different kinematical regions which may be useful for evaluation of the saturation scale. 
However, we realize that our calculations in the asymptotic regions presented in this Section are approximate 
and an exact numerical analysis of the cross section in Eq. (|39| ) has to be performed to enable us to describe 
the data with it. 



VI. CONCLUSIONS 

We have constructed an expression for single inclusive gluon production cross section in DIS (Eq. (|30|)) 
including multiple pomeron exchanges in the form of nonlinear evolution (Eq. (^)). The cross section may be 
used to describe (mini)jet production in DIS. For the case of a large target nucleus the resulting production 
cross section can be written in /c^^factorized form of Eq. ( |39| ) . The transverse momentum spectrum given 
by the cross section of Eq. (pO) reproduces the usual perturbative behavior ~ l/fc"* in the large k± limit (see 
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Eqs. ( p5|) and (|57|)). In the small transverse momentum region the spectrum softens to ~ 1/fc (see Eqs. 
( p8| ) and (^9|)). The cross section still exhibits some residual infrared divergence. This is due to the fact 
that we are scattering a point-like probe (virtual photon) on the target. The target nucleus' wave function 
has reached saturation and this is why the k± dependence softens in the infrared. However, for the cross 
section to be infrared safe (up to logarithms) , similar to |^ , we need the wave functions of both colliding 
particles to reach saturation. In pT[ | that was reached for the case of nucleus-nucleus scattering where both 
nuclear wave functions were in the saturation region. In our case the wave function of the quark-antiquark 
pair has not reached saturation yet. The onset of saturation in the 7* wave function would again be due to 
multiple pomeron exchanges. Unlike the case of a nucleus with nucleons there is no extra non-perturbative 
color charges in the original incoming 7* wave function to facilitate saturation. All the color charges have 
to be generated perturbatively. Therefore the multiple pomeron exchanges in the 7* wave function would 
take the form of pomeron loops Summation of pomeron loops is a separate problem and is beyond 

the scope of this paper. The effects of pomeron loops can be safely neglected in the energy range considered 
here. 

The result of Eq. ( |30| ) could be generahzed to the case of proton- nucleus scattering (pA). If one models 
proton as a dipole made out of diquark and a quark the inclusive gluon production cross section can be 
obtained by appropriately modifying the linear evolution term (Eq. (|l^)) and changing the 7* wave function 
into the proton's wave function. In a more generic case of a proton consisting of Nc valence quarks the 
generalization is probably a bit more involved, though the preceding evolution still remains linear and the 
gluon production mechanism would not be modified. 

We have shown that the evolution in the nuclear wave function (at least in the large-A^c limit) can be 
included via a simple substitution of the Glauber-Mueller amplitude by the dipole-evolved amplitude as 
presented in Eq. (^2|). We may therefore try to apply this principle to nucleus-nucleus collisions (AA). In 
pl| it was argued that the multiplicity of produced gluons in a central nuclear collision is likely to be given 
by the following formula: 
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(5q^]^(6) and Qos2(^) ^"^^ ^^'^ saturation scales in each of nuclei taken at the same impact parameter since 
the collisions considered are central. Inspired by the success of the substitution (22) in incorporating the 



quantum evolution corrections in DIS we may conjecture the following ansatz for the multiplicity distribution 
of the produced gluons in AA including nonlinear evolution in the wave functions of both nuclei 
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where Y is the total rapidity interval and Ngi and Ng2 are adjoint dipole forward scattering amplitudes 
in the first and second nucleus correspondingly. At the moment we can not prove Eq. (p3) and leave it as 
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an ansatz inspired by the properties of the inclusive cross sections in the saturation region studied above. 
We may argue that the final state gluon mergers which are not included in Eq. are not likely to 

give logarithms of energy and therefore should not contribute to the quantum evolution. They might be 
neglected compared to the evolution effects included in Eq. (|6^). Eq. (|63|) , together with Eq. (||), may be 
used to describe the emerging RHIC data similar to how it was done in ||41| . 
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